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Abstract. We studied the frequency resolved energy 
spectra of Cyg X-1 during the standard low (hard) spec- 
tral state using the data of the Rossi X-Ray Timing Ex- 
plorer. We found that the relative amplitude of the re- 
flection features - the iron fluorescent line at r^ 6.5 keV 
and the smeared edge above ~ 7 keV - decreases with 
the increasing frequency. In particular we found that the 
equivalent width of the iron line decreases above ~ 1 Hz 
and drops twice at frequency of ~ 10 Hz. 

An assumption that such behavior is solely due to a 
flnite light crossing time of the reflecting media, would im- 
ply the characteristic size of the reflector '^ 5 x 10® cm, 
corresponding to ^ 150i?g for a 10A/q black hole. Alter- 
natively lack of high frequency oscillations of the reflected 
component may indicate that the short time scale, r^ 50- 
100 msec, variations of the primary continuum appear in 
geometrically different, likely inner, part of the accretion 
flow and give a rise to a significantly weaker, if any, re- 
flected emission than the longer time scale events. 
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1. Introduction 

Over the last several years it has become widely accepted 
that the Galactic X-ray binaries exhibit Compton reflec- 
tion features in their spectra. 

Qualitatively, the reflection component in the X-ray 
spectra of these systems can be described as a broad 
"hump" at energies 20-30 keV, fluorescent Fe line at 
energies 6.4-6.7 keV (depending on the ionization state 
of the reflecting medium) and an absorption edge at 
the energy «7.1 keV (Basko, Sunyaev & Titarchuk, 1974, 



George & Fabian, 1991). For all geometries which do not 



obscure the direct primary continuum from the observer, 
the detected energy spectra consists of this direct com- 
ponent (approximately a power law in the 3-13 keV en- 
ergy band) and the reflection features. For standard cos- 
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mic abundances the expected equivalent width of the flu- 
orescent Fe line is ~100-200 eV (for a source above semi- 
inflnite slab of neutral matter) . 

The geometry and mutual location of the source of pri- 
mary continuum and reflecting medium should affect both 
the equivalent width of the fluorescent line and the char- 
acter of it's variability. In particular the finite size of the 
reflector implies that the time variations of the reflected 
radiation should be smeared out on the time scales corre- 
sponding to the light crossing time of the reflector. In ad- 
dition, a time lags between different emission components 
might apprear. Alternatively the geometrically different 
region of the main energy release zone (where the primary 
continuum is produced) may have different efficiency for 
production of the refiection component. As a result timing 
properties of the reprocessed component may be linked to 
the properties of the selected region of this zone. 

2. Observations and data analysis 

For our analysis we chose publicly available data of the 
Rossi X-ray Timing Explorer observations P10238 per- 
formed between Mar. 26, 1996 and Mar. 31, 1996 with a 
total exposure time ~70 ksec (we used only observations 
were all 5 PCU were turned on). For the frequency re- 
solved spectral analysis we used PGA data in the "Generic 
Binned" mode, having ji;sec~16 msec time resolution in 
64 energy channels covering the whole PGA energy band 
(B_16ms_64M_0_249). For analysis of the averaged spec- 
trum we used the "Standard Mode 2" data, having twice 
as many energy channels in the energy band of our interest 
(3-13 keV) than the "Generic Binned" mode. 

For the obtained spectra we constructed the response 
matrix using the standard tasks of FTOOLS 4.2 package 
with PGA RMF v3.5 (Jahoda 1999). The comparison of 
the timing mode spectra with that obtained from "Stan- 
dard Mode 2" showed, that they do not differ more than 
by 1-1.5%. The background was calculated using "Q6" 
model {VLE-hased model, preferable for bright sources 
does not work for these observations). However, the PGA 
background is negligible for the time averaged spectra of 
Gyg X-1 in the energy band 3-13 keV. The PGA back- 
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ground contribution to the frequency resolved spectra at 
the frequencies 0.03-30 Hz is even less important. 

3. Frequency resolved spectral analysis 

We calculated the Power Density Spectra (PDSs) in each 
energy channel for every 128 sec time segments using the 
standard Fast Fourier Transform (FFT) procedure and 
adopting the normalization of Miyamoto et al. (1991): 

P,=2\a,\'lN^R 

fe=i 
where t^ is the time label of the fcth bin, Xk is the num- 
ber of counts in this bin, iV-, is the total number of pho- 
tons and R is the average count rate in the segment 
(i? = N^/T, T - total duration of the light curve seg- 
ment). This PDS normalization ((rms/mean)^/Hz) allows 
one to obtain squared fractional rms by integrating the 
power over the frequencies of interest. Then the obtained 
power spectra were averaged over the light curve segments 
and the logarithmic frequency bins. For every obtained 
frequency bin the frequency dependent spectra were con- 
structed according to the formula: 



RXTE/PCA 



Cyg X-1 



S{E.,J,) = R,JP.if,)Af, 



2|..P^. 



T 



here S{Ei, fj) is the count rate of the spectrum on the fre- 
quency fj in the energy channel Ei. Since the light curves 
of Cyg X-1 in the different energy bands almost perfectly 
coherent and have practically the same form of the PDS 
(see e.g. Nowak et al. 1999), the spectra obtained using 
the procedure above can be used to describe the real en- 
ergy spectrum of the source X-ray flux variations on the 
different time scales. The above expession also can be eas- 
ily adapted to describe the energy spectrum of different 
PDS spectrum components. 

4. Results 

The energy spectra of Cyg X-1 in several narrow fre- 
quency ranges between 0.03 and 30 Hz are shown in Fig. |l|. 
The change of the spectral shape with the decrease of the 
characteristic time scale can be clearly seen. The relative 
amplitude of the reflection features, in particular that of 
the broad line at ~ 6-7 keV and the "smeared edge" above 
~ 7 keV, is apparently decreasing with the increase of fre- 
quency and vanishes above ~ 20 Hz. Along with becoming 
less "wiggled" the spectra become harder as frequency in- 
creases. 

In order to quantify these effects we fit frequency re- 
solved spectra in 16 frequency bins covering 0.03-32 Hz 
range in the 3-13 keV energy band with a simple two com- 
ponent model consisting of a power law and a Gaussian 
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Fig. 1. The ratio of the energy spectra of Cyg X-1 in dif- 
ferent frequency bands to a power law model with photon 
index a = 1.8. Spectra corresponding to 0.03-0.05 Hz and 
23-32 Hz were rescaled for clarity. 



line with the fixed energy En^c — 6.4 keV and the width 
FWHM=0.1 keV (that is comparable with the PCA en- 
ergy resolution at this energy). The best fit parameters 
are given in Table ^. The dependence of the line equiva- 
lent width upon frequency is shown in Fig. y. As is seen 
from the Fig. g, the equivalent width of the line starts de- 
creasing above ^ 1 Hz and falls by a factor of 2 at ~ 10 
Hz. No significant flux in the line was detected above 15 
Hz with the 2cr upper limit on the equivalent width in the 
15-32 Hz frequency range of «40 eV. 

Variations of the parameters of the spectral model, in 
particular the change of the line centroid to 6.0 and 6.7 
keV and increase of the line width to 1.0 keV, do not 
change the general trend. These variations, however, affect 
the particular values of the equivalent width and, to some 
extent, the shape of the curve in Fig. g. The change of 
the line centroid energy to the 5 keV or 9 keV, on the 
contrary, results in disappearance of the effect. 

Finally we should mention that similar behaviour of 
the frequency resolved spectra was detected in Cyg X- 
1 observations 20175-01-01-00, 20175-01-02-00 and in aU 
available observations of the set P30157 (01-10), per- 
formed in 1997 but with lower significance. 

5. Discussion 

A time averaged energy spectrum of Cyg X-1 (e.g. 
Ebisawa et al. 199^ , Gierlinski et al. 1997) as well as fre- 
quency resolved spectra at low frequency (Fig.|I|) show 
clear deviation from a single slope power law spec- 
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turcs are well known (Basko, Sunyaev & Titarchuk, 1974 



Frequency, H/ 

Fig. 2. Dependence of the equivalent width of the fluores- 
cent Fe line on the frequency. For the spectral approxi- 
mation the powerlaw-f gaussian line model was used (3-13 
keV energy band, the centroid energy and the width of 
the line were frozen at the values 6.4 keV and 0.1 keV 
respectively) . 

Table 1. Best fit parameters of the approximation of fre- 
quency dependent spectra of Cyg X-1 with the power law 
-|- gaussian line (line energy and width were frozen at the 
values 6.4 keV and 0.1 keV respectively) in the energy 
band 3-13 keV. 



^rr 



PL slope,Q EW, eV 



Freq., Hz 
(3.9 ±1.2) X 10 
(7.0 ±1.9) X 10-2 
0.12 ±0.03 
0.19 ±0.04 
0.30 ± 0.07 
0.46 ± 0.09 
0.7 ±0.14 
1.07 ±0.22 
1.63 ±0.33 
2.48 ± 0.50 
3.74 ± 0.76 
5.63 ± 1.14 
8.48 ± 1.71 
12.7 ±2.6 
19.2 ±3.8 

27.5 ±4.5 

^ - 2ct upper limit 



triim. This deviation is commonly ascribed to the re- 
flection from an optically thick media located in the 
vicinity of the production site of primary hard X-ray 
radiation. In a simple case of reflection from an opti- 
cally thick cold neutral medium the main reflection fea- 



1.93 ±0.01 
1.91 ±0.01 
1.91 ±0.01 
1.89 ±0.01 
1.89 ±0.01 
1.89 ±0.01 
1.89 ±0.01 
1.86 ±0.01 
1.83 ±0.01 
1.81 ±0.01 
1.80 ±0.01 
1.75 ±0.01 
1.70 ±0.01 
1.63 ±0.01 

1.62 ±0.01 

1.63 ±0.01 



101 ±15 

104 ± 12 

101 ±9 

101 ±8 

98.6 ±6.4 

99.1 ±5.6 

95.3 ±4.9 

87.3 ± 4.4 

78.4 ±3.9 
71.6 ±4.0 
79.6 ± 4.7 
65.1 ±5.8 
62.9 ± 7.9 

46.2 ±12.8 
<48" 
<57'" 



peorge fc Fabian, 1991 ) - a narrow unshifted iron Ka hnc 
at 6.4 keV with an equivalent width of ^ 100 eV, an ab- 
sorption edge at ~ 7.1 keV and a reflected continuum 
peaked at ~ 20 — 30 keV. However, the real spectra of 
compact X-ray binaries show considerably more compli- 
cated behavior. In particular the centroid energy of the 
line is often different from 6.4 keV, the line width in many 
cases is as large as ~ 1 keV and a broad "smeared edge" 
above ~ 7-8 keV is observed instead of a sharp absorption 
edge at « 7.1 keV ( Ebisawa et al. 1992 ). This departure 
from a simple reflection model hints at a complicated ion- 
ization state and/or motion (e.g. Keplcrian motion in the 
disk) of the reflecting media. 

It is obvious, that the simple spectral model used for 
fitting the data in the previous section is neither adequate 
nor completely justified from the physical point of view. 
The apparent line centroid energy and width vary with 
frequency. In particular the best fit centroid energy shifts 
below 6.0 keV and the width of the line exceeds 1 keV as 
frequency increases. However, due to the lack of a satisfac- 
tory realistic model of reflection in an X-ray binary, espe- 
cially applicable to the frequency resolved spectral data, 
we restricted ourselves to the simple two component model 
described in the previous section. The main purpose of us- 
ing this spectral model was to quantify the major effect 
which is clearly seen in Fig.|l|, namely the suppression of 
the reflection features in the energy spectrum with the 
increase of frequency. 

The results of the analysis presented in this paper 
show that the relative amplitude of the reflected com- 
ponent variations is lower on the short time scales (^ 
50-100 msec), than that of the primary X-ray radiation. 
Similar effect was found Oosterbroek et al. (1996) for GS 
2023±338 but the characteristic time scale in this case was 
much longer, ~200 sec. GS 2023±338 is peculiar in many 
ways. In particular the source exhibited strong (^ 10^^ 
cm"^) and highly variable on the time scales of ' 



10^ 



low energy absorption (e.g. ^ycki. Done fc Smith 1999 ) 
and, most importantly, extremely strong Fe line with 



equivalent width up to 1.4 keV ( [Oosterbroek et al. 1996 ), 
that could mean that geometry of the reprocessing media 
in GS2023±338 is different from that in Cyg X-1. Re- 
cently we found suppression of the reflected component 
variations in the low spectral state of another black hole 
candidate GX339-4 (Revnivtsev, Gilfanov & Churazov, 
1999). In this case the characteristic time scale (~50-100 
msec) was similar to that in Cyg X-1. We therefore can 
tentatively conclude that such a behavior might be a com- 
mon feature of the black hole binaries in the low spectral 
state. 

A most straightforward explanation of this effect would 
be in terms of a finite light crossing time of the reflector 
Trcfl ^ licfi/c, caused by a finite spatial extend l-^cn of the 
reflecting media. From Fig.g one can see that the equiv- 
alent width of the iron line drops by a factor of two at 



Revnivtsev M. et al.: Frequency resolved spectroscopy: Cyg X-1 



the frequency /1/2 ^ 10 Hz. This value gives us a rough 
estimate of the characteristic response time of the reflec- 
tor Trcfl ^ l/27r/i/2 ^ 15 msec and the characteristic size 
of the reflecting media /rcfl ~ 5 x 10^ cm which would 
correspond to ~ 150i?g for a IOMq black hole. We note 
however that if the primary continuum originates within 
the '^ 10i?g sphere in the inner zone of the accretion disk, 
then only a small fraction of the emitted hard radiation 
has the possibility to be reflected from the accretion disk 
regions with R <; 150i?g (in the case of a flat disk) and 
it can hardly provide observed ~ 100 eV equivalent width 
of the iron line. Therefore the assumption that high fre- 
quency variations of the reflected component are caused 
by the finite light crossing time of the reflector implies 
interesting constraints on the geometry of the source of 
primary continuum and reflector. 

An alternative explanation might be that the short 
time scale, ^ 50-100 msec, variations appear in geometri- 
cally different, likely inner, part of the accretion flow and 
give a rise to significantly weaker, if any, reflected emission 
than the longer time scale events, presumably originating 
in the outer regions. This might be caused, for instance, 
by a smaller solid angle of the reflector as seen by the short 
time scale events and/or due to screening of the reflector 
from the short time scale events by the outer parts of the 
accretion flow. 
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